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1 Introduction 

In the framework of perturbative QCD at large and low x, the conformal (global) 
group of transformations in the transverse coordinate space plays a crucial role. 
Indeed Lipatov and his collaborators [1] have derived the master equation for the 
derivative of the gluon structure function with respect to x by resuming the leading 
{cis (Q^) log(l/3;)) terms at fixed (Q^) • This derivative is a convolution in the 
transverse space of the gluon structure function times a conformal invariant kernel 
(BKKL kernel). All the relevant observables are then expressed as an expansion 
over the basis E^''^ of the conformal eigenf unctions of this kernel. The integer n is 
the conformal spin and iv corresponds to a continuous imaginary scaling dimension. 
Among others, the elastic off mass shell gluon-gluon amplitude corresponding to the 
exchange of the bare QCD hard Pomeron [1, 2] and the conformal invariant triple 
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Pomeron vertex [3]. For this purpose one is led following [1] to consider the eigen- 
functions in a mixed representation which is obtained by a Fourier transform in 2 
dimensions of E^'''". To be more specific, these eigenfunctions in the coordinate space 
are the product of a holomorphic times an antiholomorphic function. Recently, it has 
been shown [2] that the answer is a sum of two products of a holomorphic times an 
antiholomorphic form. This was accomplished by noticing that the Fourier transform 
is solution of a set of two linear differential equations one with respect to the complex 
variable and the other with respect to the complex conjugate variable. This structure 
of conformal blocks is already present in the computation of correlation functions in 
2 dimensional conformal invariant quantum field theories [4]. Various theorems have 
been established for an integrand with integrable singularities at 0, 1 and oo. In 
reference [2] a similar result seems to apply, as shown by the existence in this case 
of 2 differential equations. In the first part of this paper we show that this is not a 
sheer coincidence but that similar theorems can apply for the bidimensional Fourier 
transform. As far as the triple conformal invariant Pomeron is concerned, specific 
analytic calculations have been done recently [3]. It turns out that this quantity is a 
peculiar case of a more general calculation involving a p-tuple conformally invariant 
integral in the complex space at a particular complex value of z (resp z). The triple 
Pomeron corresponds to the case z = 1 and p = 3. Our second part is the derivation 
of this very general result. The key point is to find the set of two linear differential 
equations with respect to z and z of order p+1 which are obeyed by the p-tuple inte- 
gral. The final answer will thus be a sum of conformal blocks, as expected, each block 
being the product of one of the p + 1 solutions in z of the differential equation times 
the corresponding one in z, in order to preserve the single valuedness of the solution. 
In our study, the system of linear equations is the one obeyed by the hypergeometric 
function p+iFp{z). The case of degeneracy which is relevant for the calculation of the 
triple Pomeron vertex will be treated explicitly in a forthcoming paper. 



2 Bidimensional Fourier Transforms 



We begin this section with the following. 



Lemma 1 For < x < y < ^, let 




Then < C,R'-\ 



for any < e < 1. 
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Proof. It follows from the mean value Theorem and Jordan's lemma that for < x < 
y < |, siny — sinx > |(| — — x), consequently, 



ry 
Jo 



^ 1 _ e-W-y)y 
2R (f-y) • 



Substituting the above result into Ii, then integrating by parts yields 

TT 

h<Jj ln(| - Z/)e-^(t-^)n| - 2y)dy. 

Since |(| - |/)^ln(| - |/)| < -ft'e < oo for < y < |, e > with monotonically 
decreasing we see that 

r- 

- 2^^^4^ io ^ '^^+2^V^ (--y)-^e--(2-2/)2/d2/. 

Thus extending the region of integration in the first integral to the full semi-axis then 
doing same thing to the second integral after the change of variables yi = \ — y yields, 

which gives the result. 

With the above lemma we can now factorize certain double integrals. We will 
consider integrals of the form 



f{t)g{t)d% 

C\[0,oo) 

where C\[0, cxo) is the complex plane cut along the real axis from zero to infinity. We 
will make the following assumptions about / and g. 

i) / converges. 

ii) f{t) = f{t)e'i-\ g{t) = ~g{ty', q G »• 

iii) f{t) and g{t) are analytic for t e C\[0, oo). 
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iv) kf — f(x)dx\ < oo, kg — g{x)dx\ < oo. 



v) For e small enough but positive liniij^oo kg^R)^ = 0, linii^^oo kf{R)R'' = and 
liuiji^^ kg{R)kj={R)R^'^^ — 0, where ky{R) = maX(9g[o,27r) \g{R£''^)\ and 

kfiR) = max,e[o,2.] \fiRe'^)\. Finally 

vi) linir^o kj{r)r = = linir^o kg{r)r. 

With the above assumptions / = Imiji^^ Ir where 

Dr\{Q,R) 



with Dr being the disk of radius R. We proceed to cut -Di?\[0, oo) into the components 
indicated in Figure 1 below. If we choose the point Zi as the base point then Stoke's 
Theorem implies that 



^ JdDi+dD2 

where dDi, i = 1,2 indicate the boundary of D^, i = 1,2 respectively, = 
Jy^^ f{y)dy. Fi t being the path followed from the point zi to t along the contour 
indicated on dDi or dD2. Note that the contributions to Ir along F8541 on dD2 
cancels with the contribution of F1458 on dDi. The integral 







/■7r/2 r 


/ fi{t)g{t)dt 




/ dci> / 






Jo J 4 



f7r/2 



r-K/2 r-K/2 

< k~g{R)kj{R)R^ / d(t) j d0e 
Jo Js 



-qR{sva0—sm. (j)) 



which by Lemma 1 is bounded by 



< kg{R)kj{R)R^+'. 



(i; 



e small but positive. Since 



hit)g{t)dt 



ri8 



37r/2 

I 

7r/2 



def{Re'^)g{Re^'t')e'^^-^^R^ 



7r/2 
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An argument similar to the one above results in a bound similar to (0). The integral 



hit)g{t)dt 



/ (/7 + h-f7 + h-h + f2-h + h- f2)g(t)dt 



and is bounded by 

< 2kg{R)kj{R)R^+' + 2kfkgR' + rkj{r)kg{R)R 
The above inequalities also show that 



lim 



h{t)g{t)dt 



res 



h{t)g{t)dt- / f2g{t)dt 

732 JTi2 



0; 



and 



lim 



h{t)g{t) dt - iih - h) + {h - A) + fi{t)) g{t) dt 



This leads to, 

Lemma 2 Suppose f and g satisfy i)-vi) i/ien 



lim 

ii— >oo 



'^32 

We now consider the integral 



J«-hmg^ ((/2-/3) + (/6-/7))^?(t)rft 




c 



where C is the complex plane cut so as to make Int well defined. By rotating and 
scaling the coordinate system we may write 
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where 

C\[0,oo) 

In the above integral the branch of the logarithm selected is such that < argi < 27r 
and Int is positive for i > 1. With t — re*^ set 

and 

Thus we find 

fCXD /'2tt 



poo pzn 

Jo Jo 



Since Jq^ e''^'^ 'y)e^2ircose^Q _ q ^.i^ ggg ^j^g^^ f converges uniformly for < 

Re(7 + 7) < 1/2 and thus defines an analytic function of 7 and 7 for the values of 
these variables restricted to this domain. 
We now consider the integral 



/ / fmt)dH 



Ir = 

Dr\{0,R] 

with < Re 7 < 1/2, < Re7 < 1/2 and < Re(7 + 7) < 1/2. If = T"^ 
and gi{t) — the conditions of Lemma 2 are satisfied and taking into account the 
increment in the phase due to the cut we find 

1 poo poo 

i = Urn Ir=-(1- e'"*^) / x^-'e'^'dx / y^'^'Hy 

2 ^ Jo Jo 

where 



poo poo 
Jo Jv 



With the first two integrals replaced by their representation [6, Formula (33), p. 12] 
in terms of gamma functions we find 

/ = -(1 - e2-^T)e-2"*^e*(^+^)"/2r(7)r(7) + -(^e^^^^.^-^) _ 1)/^. 
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The first term can be analytically extended to 7, 7 7^ 0, —1, —2, .... To solve the 2nd 
term we break it up into three integrals over the following regions Ri = {{v,u).0 < 

V < l,v < u < 1}, R2 = {{u,v),0 < V < 1,1 < u < 00} and R3 = {{u,v),l < 

V < oo,v < u < 00}. It is easy to see that Ji can be analytically extended to 
7 7^ 0, -1, -2, . . ., 7 + 7 7^ 0, -1, -2, . . . and h = v^e^'^du u'-^e^'^du can be 
extended to the region Re 7 < 1, 7 7^ 0, —1, —2, For Re 7 < 1, 



/3I = 



^-^e'^ / u^-^e'^'dudv 



which converges uniformly for Re (7 + 7) < 1. Consequently I3 can be extended 
to Re 7 < 1, Re(7 + 7) < 1. This implies that / can be extended toe the region 

Re(7) < 1, Re(7 + 7) < 1, 7, 7 7^ 0, —1, —2, For Re 7 > 7(7, 7) can be written 

as 

poo poo 

7(7,7) = / / u^~^v'^~^e'^''+''Uudv 
Jo Jo 



^0 



which allows an analytic extension for Re 7 > 0, < Re(7 + 7) < 1 since this implies 
that Re7 < 1. Consequently we have shown 

Theorem 1 For < Re(7 + 7) < 1, 7, 7 7^ 0, -1, -2, . . g 7^ 

I = -(l-e^"*^)e"=^"^^g~^g^V(^+^)"/2r(7)r(7) 
+i(e2-^(7-7)_i)j(^^^). 

In particular if 7 — 7 = n, n an integer the integral becomes 
I = sin(7r7)e^('^/2)n^-7^-7r(7)r(7 + n) 

The next integral we will consider arises in the context of the QCD Domain. The 
integral to be considered is 



-II 



dh 



u+1/2 ^ V"+V2' 
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where u = —Vi + iv2, u — Vi + iv2, and Re{qp) ^ 0. With an appropriate rotation 
and scahng the above integral becomes 

1 



with 



J J (^2 _ l)«+l/2(^2 _ l)«+l/2' 

C 

and qi = qf- Since qiz comes in symmetrically with its complex conjugate we may 
assume that Reqi > and write qi = qie"^. With this notation C in the above 
integral is the complex plane cut from [l,ooe*'^) and from (— ooe~*'^, — 1]. We use 
the branch of the logarithm so that the phases oi z — 1 and z + 1 are equal to zero 
for 2; > 1. Elementary methods including stationary phase shows that the above 
integral is convergent and defines an analytic function in the variables u and u for 
— 1/4 < Rq{u + u) < 1 and Reqi ^ which is also continuous in qi in this region. We 
restrict u and u to the region —1/4 < Rert, Reii < 1/4 and consider the integral 



-// 



^2^gl/2(gi-2+qi-2) 



Dr 



_ l)«+l/2^^2 _ lY+1/2 



Using Stoke's Theorem we integrate over dDi and dD2 given in Figure 2 with / = 
(^2 _ i)-{«+i/2) and -g^iz""- 1)^+1/2 to find, 



I 



^R=of Mz)g{z)dz. 

^ JdDi+dD2 

The conditions on u and u are such that hypotheses i - iii,v, and vi of Lemma 2 
follow almost immediately. To verify iv) for the above functions we begin by showing 
that I Jp_^^ g{z)dz\ < 00. Thus we need to show that lim/j^oo I Jq" 9{z)dz\ < 00 with 

z = 1 + re**^, and with tq and determined by the equation 1 + roe*"^ = Re""^. The 
conditions on ro and (p imply that ro = R{1 + o(l)) and = + 0{^). Now 

I [ gdz\ < Cie^^^''<^-^^ I r-^^-'^dr < 00, Re^i < 1/2, 
io Jo 

and using integration by parts, 

ggi sin(0-(/i) ^-(2Reu+l) ^ 



g{z)dz\ < C2 ^ ^ j— i < 00, — - < Reu, 

|icos(0 — 0) + sin(0 — 0)1 \q\ 2 
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for constants Ci and C2 independent of R. Applying the above reasoning to / on T^^ 
shows that for — 1/4 < Re Re < 1/4, iv is satisfied. Thus Lemma 2 (appropriately 
modified) shows that 



lim - 

R—*oo 2 



J T54 J Tsi 



Ug{z)dz+ / f29{z)dz 



T32 



+ lim - 

il— >oo 2 



/ (/lo - h)9{z)dz + / f7g{z)dz + / fiog{z)d. 



Prom [7, p. 167 eq (6)] we see that 



and 



lim (/5 - h) 

K—*oo 



lim (/lo - /t) 

it— »oo 



r(i.+ |) 



Furthermore for u restricted to the region of interest the radius of small circular 
regions around ±1 may be taken to zero so that 



and 
Thus 

h 



where 



r /-oo gi5ie''*(l+r-e-*''')g-i0 

lim / q(z)dz — — ^rdr 

«-^°°Jr54 Jo (re-^'^)"+2(2 + re-*'^)"+2 

- 2ra-u) ' 



hm / W)dz = '4!^I^H^:\'li). 



7rr(l/2-M) fq^Y (11 



4r(l/2 + 'u) V2/ V2 
x[H\{q,)Hl{q,) - H\{q,)Hl{q,)\ 

+2isin7r(ii - u) sin[(7r - ^1 - qi){u - 'u)]e^^("-")/(^i), gi ^ (2) 



^''^'^ Jo ly (r(2 + re^'^))"+V2(y(2 + ye-^<^))«+V2'^'''^^- 
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Using an argument similar to the one given above but more tedious it can be shown 
that I{qi) has an analytic extension —1/4 < Keu + Keu < 1, u,u ^ 1/2 + n,n — 
0, 1, ... Thus we have shown 

Theorem 2 For u — —vi + iv2, u — vi-\- iv2, u,u 1/2 + n,n — 0,1, Ii is given 
by equation (2). In particular if vi — n/2, V2 then 

^ ^ 7rr(l/2-^) /gA» /gA" 
^ 4r(l/2 + 'u) \ 2) \2J 

[H'_M)Hl{q,) - HlM)Hl{qi)\ qi ^ 0. (3) 



3 A p-uple conformal integral and applications 

We now consider the integral 



Ip+l (ctO) <^ij 5o, Oi, bi] z,z 




We will analyze the case when a,, fij, 6, and 6, obey the following conditions: 
Condition C. For i = . . . ,p and j = 1 . . . ,p, 

a) bj, bj, Qi, di ^ Z. 

b) tti - ttj, di - dj, bijLQ - bj, bijLQ - bj, i j ^ Z. 

c) Qi — bj , di — bj ^ X 

d) Qi — di, bj — bj G Z. 
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Figure 2: 
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With = {(ii, .... cij), IP = {hi, . . . , bj), aP — (ai, . . . , CLj), V — {hi, . . . , bj), and 
Ii{ao, do, z, ^) = (1 — z)~'^°{l — z)~"'° the above integral can be recast as the iterated 
integral 

Ip+i (^ao, aF, bP; do, d^, bf; z, ^ 
= {t^ J d'zp {zpT"-' (1 - Zpf^-''^-' {z^r-' (1 - Zpf^-'^-' 

Ip(^ao, aF~^, bF~'^; do, a^~\6^~^; z Zp,z z^ , (5) 

The computation of this integral has already been done for p = and p = 1 by 
various authors [4] in the case when = di, bi = bi and gq = do- The result is 
a sum of a product of hypergeometric functions of z with its antiholomorphic part 
of argument z exhibiting the conformal block structure. Lipatov [1] has given the 
general form for p = 1 without the restrictions on the di, do, h. 

We propose to give the exact analytic structure of /p+i for any positive integer 
value of p. 

The key point of the calculation is to prove that the /p+i obey two linear differential 
equations of order p+1 namely 

OP+'{ao, oF, }f)Ip+i (ao, oF, W; do, d^. If; z, z')=0 

and 

Of+^ (ao, dP, Ip+i (ao, aP, ¥>; do, d^ , If; z,z)=Q 

where Of^^ (resp. Of^^) is the differential operator of order p+ \ defining the hyper- 
geometric function p+iFp (ao, a^, bP, z) ((resp.) p+iFp (ao, d^, bP, z^j^ 
namely 

1=1 ^ / i=0 ^ 



resp. 

1=1 ^ ^1=0 
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The general solution of this system of differential equation reads 
7,^1 {z, -z) = ^) {z) uf' (z) 



i,j=0 



where M^^^ (ao, a^, tf, z) ^resp. m^^^ ^Oq, a^, V'-, z^^ is any of the p+ 1 independent 
solutions of the differential equation 



namely 



and 



ul^^{aQ, aP, If, z) =p+iFp(ao, a*', z) , 



uf' (ao, aF, If, z) 



— z 



1-6. 



' p+iFp (oo - bj + 1, a, - + 1, 1 + bi- bj, 2 - bf, z) (7) 



Since the solution we are looking for, has to be monovalued in zz the Xij 's have 
to be diagonal provided none of the difference 6j — bj is an integer. Thus the general 
solution reduces to 



Ip+i (ao, aP, V; ao, a^, W; z, zj 
p 

= J]Af ^ (ao, a^. If; do, d^, 6^) i^f ' (ao, a^ If; z)^+' (do, d^. If; z) 

j=0 

(8) 

which exhibits a conformal block structure. 

The p + 1 unknown constants (ao, a^. If, do, d^, If^ are readily obtained by 
identifying the different behaviour of the solution /p+i near z = 0. 

At this stage we consider a generalized Euler function ([4]), 
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which is defined and analytic for Re{a+a) > 0, Re{P+P) > 0, and Re{a+a+P+P) < 
2. Using arguments similar to those given in the previous section the above integral 
can be evaluated so that for a — ci, (3 — (3 integer, 



(9) 



^ _T {a)T {j3) ^ ^^)^ sinna simiP 

- r(a + /3) T(a + i3) sin7r(a + ^)' 

Note that since a — a and (5 — $ are integer, 

sin vrct sin Tij3 _ sin ttq; sin TTf3 



sm TT a + 



siuTT (a + [3) 



To find A^+\ set 2; = in (4) in which case the integral is readily obtained since 
it factorizes. 
Thus, 



A^^ (ao, If, ao, ^) = JI a. S.-a, 

1=1 

which is ao, ao independent. 



(10) 



To find A?+\ set z. 



n 



t and observe that 



dzj (zA--' (1 - z,f'-"-' 



1- n 



,i=l 



= (di) i"^-' (1 - t) 

In a neighborhood of 2; = we find 

k=l,k^j 



W Zk \ z -t 
Kk=l,k^j 



-ao 



n (^^)' 

k=l,k^j 



-6,- 



7^+1 ~ z^-''^ (^)'-'^ ( ^ ) / (i^i ^"^-2 (1 - t)-"" (1 - i)" 



-ao 



1 _ Zk) {^_\y3-''3-'>-3+°-j 
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Thus, 



xTT-B .^1. ~7;^i7,~ (11) 

1=1, i^j 



We now give some lemmas which will help prove that Ip+i has the form indicated 
above. 
Define, 

Hjiao, If, z, z,) = {z,f--\l-z,f^-''--'u^^{ao, }f-\ zz,) , 

and 

Xj(ao, aP, If, z, Zp)^{zpY--\l-z;f--''-+'u'^{a,, aF-\ lf-\ zZp) . 
Lemma 3 Suppose condition C is satisfied. Then, 

Hj,'^ (ao, a^, F, z, z,) = (a,) ^K^^ («o + 1, a^' + 1, If, z, z,) , (12) 

1=0 P 

and for j = 1 . . .p — 1, 

0^+' H^^ (ao, aP, If, z, z,) (13) 

d_ 
dz„ 



3 

p-1 



-11 («i - + 1) («o + 1, a^ + l, V. z, Zp) 

i=0 ^ 



Proof. Since 



+ 6j - 1^ p+iFp {ai, hi, z) = {hi - 1) p+iFp (a^, hi - 1, 2;) , 
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and 

— p+iFp (oj, foj, ^) = «o n p+iFp {tti + 1, bi + 1, z) , 
^ ifc=i 
equation (12) now follows after some algebra. Equation (13) similarly follows from 
the identities, 

(^z-^ + z^~''^ p+iFp {tti -bj + 1, 1 + bi- bj, 2 - bj z) 
= {ai - bj + l)z^-''^ p+iFp {ai -bj + 2, l + k- bj, 2 - bj, z) 

+ - 1^ P+iFp («i -bj + l, 1 + bi- bj, 2 - bj, z) 
= {bi - bj)z^~''^ p+iFp {ui -bj + 1, bi-bj,2- bj z) , 
Q \ 

z— + bj-l\ z^-^^ p+iFp {{ai -bj + 1, l + h- bj, 2 - bj, z) 
(ao - b, + 1) ULM - h + l)^ ^2-., ^^^p^ _ ^. ^ 2, 6, - 6^. + 2, 3 - b, z) , 



(2-ynLi(&.-&. + i) 



and 



d 

— z^-''^ p+iFp {tti - bj + 1, bi- bj + 1, 2 - bj, z) 

= (1 - bj)z~^^ p+iFp {ai - bj + 1, bi- bj + 1,1- bj z) , 

Let Jp+i be given for < < 1 by equations (8), (10), and (11) where we assume 
that aP and V satisfy condition C. Since for < |z| < 1 the hypergeometric functions 
in/p_|_i are defined by their series representation we seek a continuation of /p+i to 
\z\ > 1. It follows from the Mellin-Barnes representation for p+iFp(ao, a^. If, z) (Luke 
[10] p. 149) that for < arg ^ < 27r 

p 

r[T(h\ n r(afc - ttj) 

p^,Fp{ao,a^,If,z) = J2^ia,)^ - 

n r(a,) .=0 n m - %) ^ ^ 

i=0 k=l 

X p+iFp (^ttj, l + Qj - bi, l + ttj- Qi, ^ (14) 
The following Lemma will be useful it what follows: 
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Lemma 4 Suppose the points Xi and yj, i — . . .p, j — 1, . . .p are distinct and none 
are equal to one. Then, 



B. 



n,k 



P 

-TT i- — Xj 

^ (1 - Xn){l - Xk){Xo - yj){yj - Xj) -pr {Xj - Vj) 

^ (1 - yj){l - xo)(yj - Xn)(yj - Xk) jjf (yi - yj) 



J 



Ui 



UO IVi=i{Uj-Zi) 

where Zk — y^ — 1 and Ui — Xi — 1. 
Proof. Bn^k is equal to 



■K, (15) 



fj- ZjUQ^Zj - Un) (Zj - Uk) IJ- (Zi - Zj) 



Write III ft = Sn,k where, 



'3> 



,=1 zA^i - '^n){Zj - Uk) {Zi - Zj) ' 

and A = ^ nf-i -■ The function 

UQ ^ 1.1—1 Ui 

p 



^(^ - '"n)(2; - Uk) {Z - Zi) 

behaves as z^"^ for large z, and by hypothesis has simple poles at z = Zj and z — 0. 
If n = k then there is one extra simple pole at 2; = -Ufe. Cauchy's theorem says the 
sum of the residues is zero and the evaluation of the residues gives the result. 
With s{ai) = sinTraj we now prove, 

Lemma 5 Suppose ao, Oo; oP, If and If satisfy condition C then 
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ip+i{z, z) = Jjje'^^"' "^z ^'p+iFp (ttj, ttj - bi + 1,1 + ttj - ai^j, 
i=o ^ 



X p+iFp ( cij, cij — bi + 1,1 + dj — ai^j, - 



where Jjj is given by, 

s(bj - dj) 



i=0 



and 

p 

JqO ■^bi—ai,ai—ao,bi—ai,ai—ao- i^'^) 

i=l 

Proof. From equation (14) it follows that 

P+iFpiao, a^, If, z) p^,Fp{d,, d^, lf;z) = \{ ~ 

1=1 

I •J' ^ / 6*"^ \ / 



where 



and 



i=l 



c^j- = B{aj, ao - aj) ]^ B{bi - ai, ai - aj). 



i=l 



with B{a,b) = ^f^^^- c^^ is obtained from cJJ^ by replacing ao, aP and 6^ by ao, a^ 
and 6^ respectively. Likewise we find 
p 

J2 yj^'z'-'^z^-~'^p+iFp{ai - bj + 1,1 + bi- bj, 2 - bj, z) 
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X p+iFp (^ai -bj + l,l + bi-bj,2~ bj] ^ 

J^i s{bj - ao) 

Yj s{(i, - hj + l)s{h, - hj) 

A. 



n 



,=1 s{bi - bj + 1) 



X 



X 



c^fc \ V) 1 + afe - 6j, 1 + Ofe - a^^fc; - J 

^ c^fc f j 1 + Ofc - 6il + afc - ai^k, - 

i~r\ \ ^ / ^ 



with 



and 



r(2-6^r(V-i) 
^o = r(6,-a,)r(a,-., + i)n^(^-«^'«-«o), 



^ r(6j- - ao)r(ao - 6j + 1) f^f 



With the simphfications 

s{di — bj + 1) s{di — bj] 



= 1, 



s(bi - bj + 1) s{bi - bj) ' 

jib,-m2-b,)Tib,-l) 
s{bj - ao)r(6j - ao)r(ao - 6j + 1) 

s{l - dp) _ s{do) 
.(6,-1)- 

and 

s(a,/) 

5(1 - ao, ao - ai) = B{ai, gq - ai)— — r, 

s[ao) 
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the coefficient J„fe, n 7^ 7^ A; of 



p+iFp {an, an -bi + 1,1 + an - ai^n fraclz) 



Xp+iFp I Ofc, Ofe — 6i + 1, 1 + ajfc — ai^k', 



becomes, 



Jnk = Yi^i^i~^i^^i~^ri)YlB{bi-ai,ai-ak) 

i=0 i=0 
P 

X ]^sin7r(6i - ai)B{an,ao - an)B{dk,do - dk)Jnk, 



wfiere 



Jnk — Y\. 



i=l 



s[ai 



y^^ ^y,--a, s{bj - ao)s{an)s{hk)s{bj - a^) 



X 



p=i s{bi - bj) 



Now from tlie fact tliat Oj — Oj and bi — bi differ by integers we find, 



[-ly 



s{bj)s{ao)s{bj - an) 
s{bj - do)s{do)s{dn)s{dj - bj) 
s{bj)s{do)s{bj - dn) 



so tliat 



J, 



nk 



1 

^ s{bj - do)s{dn)s{dk)s{dj - bj) -A- s{di - bj) 



E 



n 



s(6j)s(ao)s(a„ - 6j)s(ajk - bj) t=f - bj) 
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For n — k the coefficient Jofe of 



p+iFp ( 1 + Oo ~ 1 + Oo — Oj, - 



X p+iFp I Ofc, 1 + Ofe — 6i, 1 + Ofe — Oi^fe, - 

z 



becomes, with the substitutions above, 



Jo,k = - ai,ai - ao)Y\_B{bi - ai,ai - ak)B{ak, Oq - ak)Jo,k 



where Jo,fc is given by (18) above with n — k. For A; = we find 

p p p 

Jo,o = Y\_B{bi- ai,ai- ao)Y\_B{bi- di,ai- di)Y\_s{bi- di) 



1=1 



X 



i=l 



s{di) s(ao)s(6j — ttj) -A- s(aj — bj) 



n 



1 ^ ^{bj)s{bj - dj) s{bi - bj) 



(21) 



With the substitutions yj — exp —2mbj, Ui — exp —2mai and = -^n.^nr expi^b' 
the result now foUows from Lemma 3 and the symmetry between J„ o and Jo,n- 
Define 



i=0 



and -0^ as above with bi and replaced by bi and respectively. 

Lemma 6 Suppose Condition C is satisfied, V'p-i is not an integer, Re('0p-i+'?/'p-i) > 
p — 1, Rc(ap + dp) < 0, Re{bp + bp — ap + dp) < 0, Re{bj + bj — Up — dp) < 2, and 
Re(6p + bp — Qj — dj) < 2. Then for z or 1 



/p-i 
dzpUo Y\ CLi Kl{aQ + 1, + 1, If, z, Zp)H^{do, a^, z, Zp) 
i=0 



(22) 



p-i 



+ lli^i - bj + 1) K^{ao + l,aP + 1, F, z, Zp)H^{do, d^, F, z, Zp) = 0, 



i=0 



where F = [0, oo exp — i0) U [1, oo exp —i^), (f) — arg(^) 
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Proof. We use F a cuts for 2;" and (1 — z)" choosing the determination for both 

so that both give positive reals when their arguments a large positive real numbers. 
We beginning by breaking up the above integral into pieces for which \zZp\ < 1 and 
\zzp\ > 1. For \zzp\ > 1 we use the fact that 

p-i 

H aiXl\ao, aF-\ lf-\ ao, a^"', (23) 

i=0 

1=1 

and 

p-i 

l[{ai - bj + l)Af '(ao, aP-\ ¥'-\ So, dP-\ If-') (24) 

i=0 

p-1 

= ao n(^i -^i- l)Ar'(ao + 1, a^"' + 1, ao + 1, a^'' + 1, 6^"'), 

i=l 

and Lemma 5 to represent the integrand in (22) for \zZp\ > 1. The conditions on 
Oq . . . , Op, 61 . . . , flp, o-o . . . , flp, and bi . . . ,bp imply that the above integral is convergent 
for Zp ~ 00,0, and 1. Furthermore with condition C and the condition on ipp-i the 
hypergeometric functions in the integrand behave as pFp_i(ao, a^, if, 2;) ~ (1 — z)'^^-'^ 
([8]) so that the integral is convergent for zzp ~ 1. The last two parts of condition 
C insure that the contributions along the two sides of each cut cancel out giving the 
result . 

We now prove our main result, 

Theorem 3 Suppose that condition C is satisfied, ipp-i is not an integer, Re('0p_i + 
ipp-i) > —2, Re{ap + ap) > 0, Re{bp + bp ~ ap — ap) > 0, Re{bj + bj — ap — hp) < 2, and 
Ke{bp + bp — aj — aj) < 2. then /p+i is given by equation (8) with A^"*"' given by (9) 
and (10). /p+i may be extended using the above representation so that only condition 
C is satisfied by the parameters. 

Proof. Consider Jp+i for p = 1, the result is already known ([4]) and we will recover 
it. In this case, 

h = (25) 
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c\r 




c\r 




^z^-'il - z^f^-''^-'zr^-\l - -z,t-~^^-'h{a^ra^,zz,,zzi) 

^Xc\r^r-'(1 - z^f--'^--^z^--\\ - z^f-~^--^ 

c 

xIi{ao,ao,zzi,zzi). (26) 

where F = [0, cxoexp — #) U [1, ooexp — (j) = arg(2) and Xc\r is the characteristic 
function of the set C\r. We will use the determination given in the above lemma 
for z". We suppose that 2; 7^ or 1. If Re(ai + Si) > 0, Re(6i + 61 — ai — Oi) > 0, 
Re(ao + do) < 2, and Re(6i + 61 — Oq — ao) < 2 the above integral converges uniformly 
on compact subsets of the z plane that exclude zero and one and defines an analytic 
function of the a and b variables. Utilizing the third line in the above equation it is 
not difficult to see that the integral is continuous in z. For e small but positive let /| 
be given by the integral above with F replaced by F^ which is the contour with every 
point in F displaced by a distance e. Thus 

/| = / / d''z,xc\TA'~\'^ - z,)'^-'''-'zr'-\l - z,)''-^'-'h{ao, do, zzi, zz^). 



c 

The integral converges uniformly for z on compact subsets contained in or on F^ that 
avoid z = 0,1 and defines an analytic function in the variables Oq, ai, bi, do, di, bi in 
the region designated above. Furthermore it follows that /| is a continuous function 
of e for z sufficiently close to F\{0, 1} and hmg^o^l — -^2 where the convergence 
is uniform on compact subsets of F\{0, 1}. Since for every e small but positive the 
integral, and the integral with the integrand differentiated once, and the integral with 
the integrand differentiate twice with respect to z converge uniformly, it is possible 
to interchange integration and differentiation with respect to z. Prom (12) above we 
find 

If 

OI 1\ {z, z) = -«o— / ^^^^'q^^o^ ("^0 + 1, ai + 1, bi, z, zi) 

H^^^ (ao, ai, bi, z, zi^ 
Thus from Stokes' Theorem we find 



Om = ^J^ K^'\z,zi)H{z,zi)dzi 
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= Ge{z,z) 



To continue on we impose the extra restriction that Re(ao + ao) < 0. Under these con- 
ditions it follows using Stokes' Theorem that for z sufficiently close to r\{0, 1}, G^{z, z) 
is continuous in e. Prom Lemma 6 we also see that limg^o Ge{z,z) — uniformly on 
compact subsets of r\{0, 1}. Thus we can conclude that for z inside T^. 

I^{z) = i/Jh{z)+M^) 

where iph{z) — CqUq + c\u\ and 

where w{ul,u1) is the Wronskian of and uf while Wi{ul,ul) is the determinant 
obtained from w{ul, uf) by replacing the ith. column of w by d") . Condition C insures 
that Uq and u1 are linearly independent. Here Zq is chosen so that < |2;o| < 1 and 
zq lies on F. The coefficients Cq and c{ are functions of z (thought of as independent 
of z for their computation). A similar discussion interchanging z,ao,ai, and bi with 
z, So, oi, and bi respectively shows that 

1 

-^2= XI (^lj'^i((^o,ai,bi,z)u'^j(do,ai,bi,z) + f, 

i=0,j=0 

where is a function of ^(oq, ai,bi,z), ^(oq, Si, &i, z) and their first derivatives. 
The theory of differential equations [9 Theorem 4.1] says that the and p are 
continuous functions of e for z G r\{0, 1} and lime_^.o/^ = uniformly on compact 
subsets of r\{0, 1}. Writing 

—^z1'-\l - Zif '-^-'-^l - zziy^zi^'-^l - ZiY'-^^-\l - 
2z 




c\r 
Nil<i 



c\r 
zi |>i 

= Ii+Il (27) 
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With the above constraints on the parameters we see that in a neighborhood of 2; = 0, 
II ^ K where A' is a constant independent of z and z. The change of variables t — zzi 
in the second integral yields 

j2 _ ^1-61^1-61 ^_-|^'^6i-6i-ai+ai f f ^ ^ 'jbi-ai-l ^-^ ^■^-ao^i-2^-[^ ^-^ 

J J 2,i t t 

C\[0oo) 
l4l>l 

Consequently the dominated convergence theorem shows that 

lim/i-if^i-^/2 = A^. 

Since the only constraints on bi and 61 are those give above and condition C the 
above considerations imply that /3q^ = = /3°q and = We also find from the 
dominated convergence theorem that for z real lim2^o-^2 = Aq. Prom the previous 
discussion the solution reads 



and 



I2 (ao, oi, 61, oo, Si, 61, z, z^ 

A[,^^ (ai, hi] ai, 6ij 2^1 (ao, ai, hi] z) J^i (ho, ai, hi, z^ 

+Xf^ (ao, ai, hi, Oq, Oi, bi^ {zf~''' {z)^~'^' 
2F1 (oo -61 + 1, oi - 61 + 1, 2 - 61; z) 
2F1 (ao -61 + 1, fii - 61 + 1, 2 - 61; z^ 

Aq^^ (ai, hi, di, bi^ 

r (ai) r {hi - ai) ^ (^1) ^ (^1 ~ ^1) ^^"^ ^ (^1 ~ ^1) 

^(^0 sin7r6i 

A^^ (ao, ai, bi, do, Oi, 

r (bi - 1) r (-ao + 1) r (^i - 1) r {-do + 1) 



= (-1) 



bi—bi+ai 



r (-ao + 61) 



r (-So + 61) 



siuTrfci sin7r(— Oo) 
siuTT (61 — Oo) 
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in agreement with previous results. This proves the result for p = 1 for Re(ai+ai) > 0, 

Rc(6i + &i — fli — Si) > 0, Re(ao + < 0, and Re(6i + bi — aQ — Sq) < 2. We may 
now extend I2 using the above representation so that only condition C is satisfied by 
the parameters. 

Suppose now that, 



then 



p-i 

= J]Af (ao, a''-\ }f-\ ao, ~oP-\ y'"^) (ao, a^-\ lf-\ z) x 

j=0 

X u^. (ao, a^-' F-\ z) , 

Ip+i (ao, aP, V, ao, oF, W, z, z^ 

p-i 

= E ^"0 ^ 

xi /" /" d2zi(zi)"^"'(l-^i)'^""^"'<(ao, a,, 6,, 



c\r 



{z,f^-\l-z,f--~^--'u^^(ao, aP-\ ~V>-\ zz,) . 

The conditions imposed above imply that Jp+i is analytic in the parameters ao, a^. If, 
clq, aP, and IP in the region specified. Furthermore Jp+i is a continuous function of z 
for 2; 7^ 0, 1. We now impose the condition that Re('^p_i + 'ipp-i) > p — 1. Let be 
as above and set 

r<: — 
-'p+i ~ 

xIp{ao, a^~\ 6^~^ao, a^"\ ^^i), 

and a discussion similar to that given above shows that limg^o Ip+i — Ip+i uniformly 
on compact subsets of F\{0, 1}. Since Ip^i can be differentiated twice we find using 
Stokes' Theorem that 
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1 zl' (1 - z^f^-"^ 



O^z^^ (^i' ^i-i ^0; So, ^, z) 

P-1 ^ 

-do n ~ ~ 2i / / '^^^ 

c\r, 

Ip [aQ + 1, + 1, 6i, oo, Oj, 6i, zzi, zz^ dzi zp''^ (1 - zif 



+ K;(ao + 1, a^' + 1, fe^', z, z,)H^{ao, a^, z, Zi) 

Lemma 5 shows that the last integral is convergent and lim^^oG^p+i = uniformly 
on compact subsets of r\{0, 1} by Lemma 6. A similar equation is obtained when 
using Of'*'^ and we are led to the solution 

p+i 

i,j=0 

where condition C insures that the Uj~^^,j = 0, . . . ,p are linearly independent. We 
again appeal to the theory of differential equations to show that limg^o P — ^ a-nd 
that is a continuous function of e. Using the above continuity properties we find 

that 7p+i = ESo Ajt^r'<^'- Write J^+i = E?=o^p+i> where, 

lU ^YiJ I d'z,x\..,\^,{z,)zr'zr\l - ^i)^— ^(1 - ^i)^-^-^ 
c\r 

xii^(ao, a^-S zzi)u^^{ao, ~aF-\ zz,), 

for J = 0, . . . , p — 1, and 

TP _ 

-'p+i ~ 

J=0 c\r 

|zzi>l 

Xp-fp-i(fli, % — 6i + 1, 1 + — a^^j, )pFp_i(a,-, a,- — 6^ + 1, 1 + a,- — Oj-^,-, 3^)- 
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Prom the dominated convergence theorem we see that hm2;_^.o Ip_^.i — Aq^^ and 
hm^-^o z^^~^z^^~^Ppj^-^ — \^^^ for z real. Furthermore with the change of variables t — 
zz\ it is not difficult to see that lim^^o z'^p~^z'^p~^Ip_^i = K^~^^, with Kp~^^ independent 
of z. Condition C and the above discussion shows that (3ij = 0,i ^ j and f3i^i ~ 
Af"*"^,? = 0, — 1. In order to compute (3p^p split /p+i up into the pieces Ip_^.l 
for \zi\ < 1 and Ip_^_^ for \zi\ > 1. With the change of variables t — zzi the second 
integral becomes, 

C\[0,oo) 

xipiao, aF-\ lf-\ fio, aP-\ t, t). 
Using the dominated convergence theorem we see that, 



^_l\bp-bp-ap+ap z'^''~^ Z^^~^ I 




dHt''^-H^^-^Ip{aQ, aP-\ lf-\ ao, aP-\ t,t). 

C\[0,oo) 

We now show that the above integral is equal to (^—i^p-bp-ap+ap xp+^ _ Note that 
the result has been shown for p = 1. Now suppose the result is true up to p. With 
the substitution of (5) for Ip we find, 

^_l\bp-bp-ap+ap ^-^ 

^2^^6j,-2^p-2 J J ^2^^ap_i-l^aj,_i-1^2 - - ^fp-i-^p-i-l 

c\[o,oo) c\r 

7p_i(ao, a^'-^ So, a^-^ tz, tz). 

Decompose the second integral into the regions Ri = \z\ < 1 and i?2 = |^| > 1 then 
make the change of variables y — zt. The integral Ki over Ri becomes, 

2^^b„-276„-2 f f ,2 fy\ fW^-'-WyV"- 





C\[0,oo) C\[0,oo) 

I |j Ip-i{ao,aP ,Jf ,ao,aP ,If ,y,y), 
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where B the closed unit ball in complex plane. Interchanging the order of integration 
then setting y — ut yields, 



C\[0,oo) 

J y"d2^(M^)''^-2(M^^^p-%_l(ao,a^'-^^^-^ao,a^'-^6^'-^^^^,M^), (28) 

c\r 

where r is the cut described at the beginning of the Theorem with = argi. Finally 
with the change of variable w — ut we arrive at, 



C\[0,oo) 
|u|<l 




c\r 



Applying the same operations to the integral over region R2 and utilizing the fact 
that for III > 1, (1 - D&p-i-^p-i-i = |6p-i-ap-i-i^^ _ ;^^6p_i-ap_i-i ^Y^en adding the 

two integrals gives, 

(_l)&p-&p-ap+Sp lijn ^bp-i#p-i/^ (30) 




c\r' 

p-2 

B ~ ~ Y[ B 

bp — I, — ao+1; bp — 1, — ag+l \_\_ a^ — 6p+l, bi—ai] at—bp+l, bi—ai'> 



i=l 



where the induction hypothesis has been. The claim now follows. 

We may now extend the result to the case when Re(-?/'p_i + V'p-i) > ~2. The 
above representation may be used to define Jp+i when only condition C is imposed 
on the parameters which completes the proof. 

A comment is in order. The result of this integral for jz] < 1 is easier to get by 
using a new basis namely the p^iGp{ai;bi;z), which are nothing but the hypergeometric 
functions p+iFp{ai;bi;z) up to a multiplicative factor namely 

p+iGp{ai;bi;z) = j-j-p r{b ) p+^'^p^'^*'^''^-' 
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The basis of the expansion is now Vj{z)yj{z) rather than Uj{z)Uj{z) where 
Vj{z) = z^-''^ p+iGp{ai - bj + l;bk - bj + l;z) 

and 



Vj{z) = z^ ''^p+iGp{di - bj + l;bk - bj + l;z) 

In these formulas k ^ j = 0, p where by convention 60 = ^'o = 1- 
In this basis the solution reads 

i=o 

where the fij are obtained from the corresponding values of the Xj. 

It turns out that an obvious factorization does occur. Namely, by defining 



1 P 

we get 



1=0 



Pp+1 T-rp 



Ulo^Sibj-k) 
To be more specific 



's{hmi^,^^s{b,-h) 



Similarly, for \z\ > 1 it is useful to expand the integral as a combination of 
hypergeometric functions which are defined for \z\ > 1 namely 



(z) p+iGp{aj,aj -bi + l;aj - + 1;-) 

z 



or more specifically 



Wj{z) = {z) p+iGp{aj,aj - k + l]aj - + 1;-) 

z 
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the solution reads now 

I^Y.v,W,{z)W,{z) 

where the Vj have to be determined from the value of the integral for 2; — > oo. 

We proceed as before and it turns out that the same factorisation does occur 
namely 

Comparing iij and Uj amounts but for a sign to interchange hk and as expected 
with the interchange of 2; — > K Indeed if in the integrand we change -Zj — > ^ the new 
integrand J reads 

p p 



J = {-iy-'z-''°z-^°Y[ti''°-'''{i-ti)''*-''*-\i-Y[tiZ-^) 

i=l 1=1 
P _ P 



1=1 i=i 

It is now obvious that the result of the integral for > 1 is deduced from the result 
for |2;| < 1 by a simple change of the parameters of the hypergeometric functions of 
argument ^ up to the multiplicative factor 

namely if we define ctj and Pi as the parameters we get by simple identification 

«i = ao-bi + 1 
(5i = ao - aj + 1 

which corresponds to the parameters of ^0(2;) and 

ai- I3j + 1 = aj -bi + 1 

2 — Pj = ttj — oq + 1 

which corresponds to the parameters of Wj{z). The above argument can be made 
rigorous using (14). 
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To summarize 

p 

Ip+i = iip+i^iijVj{z)Vj{z) 

j=0 

where 

Vj{z) = z^-''^ p+iGp{ai - bj + l;bk- bj + 1; z) 



1 P 

and 



1=0 



and a similar formula for 1 2; I > 1. 
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Appendix A 



In this appendix we give an alternative method to calculate (30). Besides the condi- 
tions imposed in the hypothesis of Theorem (3) it will assumed that bp—bp—bi+bi > 0, 
i = 0, . . . ,p — 1. We need to calculate 



i=Q 

where 

Ui{t) = t'-'''pFp^i{aj -bi + l;bk-bi + l;t) 

where j, k i = 0, p — 1, and Af~^ is given in formula 11 and by convention bo = 1 

We split the domain of integration in two parts, and write I = Ii + I2 where Ji is 
the contribution for \t\ < 1 and I2 for \t\ > 1. In the first domain the hypergeometric 
functions arc convergent so that 



1=0 



B/r 



where B is the unit disk in the complex plane and F is the contour used in the proof 
of Theorem 3. To calculate I2 , the other form of the solution which is valid for > 1 
will be used, namely 

/2 = j2''r' J I t'^~^Wi{t)t'--^Wi{t}dH 

where 

Wi{t) = r"^ pFp_^{ai - bj + l;ai-ak + 1; -) 

V 

with the same convention for the indices j and k 
By definition 

^p-i iTjzUnbj -bi+ ml, - h + 1) ^ 

' m:;r(%-6, + i)r(a,-6, + i) 



llLo - bj + l)r(ai - bj + 1) 
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p—1 p— 1 

1=1 1=0 

p-i 



TT 

1=0 

and 



For ■{/; real and \w\ <r <1 set, 



= yy {wtf^-^Ui{wt){wtf^-^Ui{wt)dH 

B/T 

and expand (resp UiiJ)) as a series (resp ^^). The integration on the 

azimuthal angle yields at once 

m = n + hp - hp - hi + hi. 

The integral over the modulus of t is also straightforward and yields an additional 
factor 

1 T{n + hp-hi) . 2^n+6,-6, 

2V{n + hp-h + iy ' 
Collecting all the pieces we get 

where 2p+iG2p{w'^) is the unrenormalized hyper geometric function the upper param- 
eters of which are aj — hi + l,aj — hi + 1 + bp — bp,hp — hi and the lower ones read 
bk - hi + l,bk - hi + 1 + bp - bp, bp - hi + 1,1 + bp -bp - hi + h^. 
Now define 



'""^ > '^''2inJc, r{l - a, -t)r{hp- hp + h,+t) r{t + bp) 
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where Cr is a loop beginning and ending at +00 and encircling all poles of r(... — t) 
once in the negative direction but none of the poles of r(... + 1). This yields at once 

hiw) = fniw) 

The same method holds for l2{w), w real and \w\ > r > 1 with the integration 
over the azimuthal angle yielding n = m + bp — bp — ai + di while the integration on 
the modulus of t gives 

i rfe-6,+ l + m) 

'2r{di-bp + 2 + my ^ 

As before hiw) is a comphcated sum of G functions. 
Now define, 



p-l 



r(l - b, - t)T{bp -b, + ~a, + t) r(t + 6,-1) 



M r(i - % - 



t)V{bp -bp + bj + t) T{t 



where Cl is a loop beginning and ending at —00 and encircling all poles of r(... + 
once in the positive direction but none of the poles of r(... — t). 
It is easy to check that 



= —hiw) + residus t — 1 — bp. 

Indeed 

f^-'^i-iy-' 

But Cauchy theorem implies that 



IVjZi^iSidj - di) 



hiw) = fniw), 
so that in the limit w — > 1 we find, 



^2 yi - bj)ril -bp + dj) 

'~'f}oribp-aj)rii-bp + bj) 



36 



or equivalently 

j_ T-i^ >S'(6p - aj) YT r(l + aj - hp)V{l -bp + dj) 

- n s{b, - h,) \\ r(i -h, + h,)T{i + \ - h,) 

At this stage it is worth noticing that 

jip = iJ,p-iT{bp - ap)T{bp - dp)S{bp - ap), 
which can be recast as, 

V{l + ap-hp)V{l + ~ap-hp)S{hp-ap) 

this yields at once 

This is the expected result. 

As a nice application we may recover the generalized Euler formula, namely cal- 
culate 

j t''-'Uo{t)f^-'Uo{t)dH 

where Uo{t) = {1 — t)^^^ =i Fo{b ~ l,t). If we follow the method presented above we 
are led to calculate 

1 f T{-t)T{a-a+l-b + t) T{t + a) 

JR = 7T^ / dt- 



2m Jc^ V{b-t)V{a-d + l + t) V{t + a + l) 

since by identification 6i = a — 1, oq = 1 — fe, it is easy to get that 

fn = IS{b)V{l - b)V{l - ly 

by picking the pole at t = —a, we get at once 

^ V{a)T{l-~a-b) 
r(6 + a)r(l-a) 

collecting all the pieces the final result reads 

T{a)T{b)T{l -d-b) 



r(i -6)r(6 + a)r(l - a) 
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